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We apply wave packet methods to study an ion-trap system in the strong excitation regime 
imposing neither the rotating wave nor the Lamb-Dicke approximations. By this approach we show 
the existence of states with restricted phase space evolution, as a genuine consequence of quantum 
interference between wave packet fractions. A particular instance of such a state oscillates between 
maximal entanglement and pure disentanglement between the constitute subsystems, where the 
characteristic crossover time is very rapid. Over longer time periods the dynamics of these states 
exhibits collapse-revival patterns with well resolved fractional revivals in autocorrelation, inversion 
and entanglement. 



PACS numbers: 42.50.Dv, 32.80.Xx, 33. 80. Be 

I. INTRODUCTION 

Progress in research fields as laser cooling, produc- 
tion and controlling sub femtosecond laser pulses, man- 
ufacturing of solid state devices, and so forth, ren- 
der experiments of novel quantum characters. A key 
model in quantum optics is the harmonically trapped 
ion pumped by external lasers [lj] , which successfully has 
been used to study pure quantum phenomena in experi- 
ments Q, S 0, . Both the internal electronic structure 
of the ion and its center-of-mass motion are treated quan- 
tum mechanically and it is possible to tune experimental 
parameters such that a single or few electronic transi- 
tions in the ion can be isolated and coherently coupled 
to the motion. The spatial profile of the driving laser, 
either a standing wave (SW) or a traveling wave (TW), 
effectively reshape the trapping potentials. In particular, 
a unitary transformation of the regular one-dimensional 
ion Hamiltonian with TW pumping may result in two 
equally shifted and displaced coupled harmonic oscilla- 
tors [y, 0] , whereas S W fields in general generate coupled 
potentials with a more complex shape. In addition, a 
fairly new proposal, verified by experiments, is to place 
the trap inside a high-Q micro cavity such that it is the 
cavity field driving the trapped ion [8|, |9|]. Even Bose- 
Einstein condensates have been successfully trapped and 
coherently coupled to a single cavity mode [10[. 

In most theoretical work, pertaining to both of the 
SW and TW cases, approximations are imposed in order 
to obtain analytical or semi-analytical results valid in 
different parameter regimes: 

1. The Lamb-Dicke (LD) regime. In this regime the 
wavelength of the classical laser field is long com- 
pared to the extent of the confining harmonic trap 
and an expansion in the small Lamb-Dicke pa- 
rameter r\ of the mode profile function is made 
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2. The rotating wave approximation (RWA) regime. 
Here, a basis is defined in which the involved time 
scales differ considerably and, once in the inter- 
action picture, fast oscillating terms are neglected 
Ell EI El US This assumes that the par- 
ticular system can be separated into a "slow" and 
a "fast" part, which puts constraints on the con- 
stituent frequencies; trap frequency lj, ion transi- 
tion frequency laser frequency lol and Rabi fre- 
quency A (laser- ion coupling amplitude). 

3. Strong and weak excitation regimes. These are the 
limiting situations of either a very large or small di- 
mensionless parameter X/lj [ll, [23|, [24[ correspond- 
ing to, respectively, strong or weak pumping of the 
ion. It should be noted that this definition assumes 
a moderate number of phonon excitations, other- 
vice the effective Rabi frequency may become rel- 
atively large even in the weak excitation regime. 
This in fact, is the case of the current work. 

Typically, the parameter regimes for which the above 
approximations are justified overlap and, in particular, 
there are regions in which none of them is valid. In 
many studies, one of the approximations is applied and 
one considers dynamics "beyond" the others. This, how- 
ever, can be misleading since one is often in a parame- 
ter regime where also the non-imposed approximations 
could in principle have been implemented. Recently, Liu 
et al. [25[ showed that the ground state of a TW pumped 
trapped ion may be considerably lowered when both the 
RWA and the LD approximation are invalid. In the cur- 
rent work we use a fully numerical method and therefore 
no approximations are imposed. We discuss, neverthe- 
less, the various parameter regimes as well as introduce 
an additional approximation, namely, the adiabatic ap- 
proximation, to provide further insight into the quantum 
dynamics. It is commonly thought that the dynamics, 
in contrast to the simple cases of the validity regimes 
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discussed above, becomes irregular beyond these approxi- 
mations. We will here show that this is indeed not always 
true. 

Theoretical and experimental research on trapped ions 
has been concerned with both SW and TW driving, 
with the main focus on the latter. Examples include 
state preparation of non-classical vibrational states such 
as Schrodinger cat states or Fock states for both TW 
[H [H EH, 13, [H and SW [lip, HJ [IH pumping, 
collapse-revivals in the TW [fll5| and SW cases, 
state measurement [l6| and quantum information pro- 
cessing d[27|. 

The above-mentioned theoretical works almost exclu- 
sively employed the formalism of vibrational creation and 
annihilation operators for the harmonic trap. We, on 
the other hand, will reformulate the model in terms of 
ionic center-of-mass position and momentum, casting the 
problem into one of two coupled harmonic oscillators. In 
particular, we will study the coupled dynamics of an ion 
wave packet evolving on the two potentials in the case of 
SW pumping. For a one-dimensional system, as the one 
we consider, it resembles an idealized diatomic molecule 
for which only two coupled bound electronic states are 
taken into account. A similar analogy, but between a di- 
atomic molecule and a cavity QED model, was pointed 
out in [H, [29[, works that also made use of the wave 
packet method. As the standard viewpoints and meth- 
ods applied in quantum optics and molecular physics are 
rather different, we discuss them and their relation in 
some depth with the aim to facilitate branching out and 
combining the two research fields. 

Since the pioneering work in the seventies by Heller 
[30( | . wave packet techniques have been used especially 
in molecular and chemical physics research [3l|. It is 
not always possible to separate the dynamics in isolated 
electronic states and the evolution on coupled electronic 
potentials has to be considered; couplings could be of e.g. 
vibronic, spin-orbit, or rotational nature or induced by an 
external laser field. In Fig. HJ we display two schematic 
examples of coupled potentials, bound-bound in (a) and 
bound-repulsive in (b). Solid potential curves are dia- 
batic, while the corresponding adiabatic ones are dashed 
(for a proper definition, see Sec. Ill Bp . When the wave 
packet propagates on the coupled potentials, interference 
effects will occur. In the bound-repulsive system, one 
such effect manifests itself as unexpectedly long-lived vi- 
brational states [32[. This takes place when the interfer- 
ence restricts the wave packet evolution in such a way 
that the repulsive part of the system at the right side 
of the curve crossing is not reached, indicated with the 
heav y so lid lines in Fig. [TJb. We call this a bistable mo- 
tion [33|, [34|, [H, [3^ . Similar interference effects can also 
occur in a system consisting of two coupled bound states, 
in which case the wave packet always returns to the same 
diabatic potential after one oscillation (see Fig[TJa). In 
astable motion (not shown), the wave packet starts in one 
diabatic state, splits at the curve crossing and when it re- 
turns to the curve crossing it recombines in such a way 



that it switches completely to the other diabatic state. 
Consequently, astable motion cannot exist in a system 
comprising coupled bound-repulsive states. We have pre- 
viously carried out a numerical wave packet study of the 
dynamics of bistable trajectories [36[ and found that the 
trajectories can be exceedingly long-lived with sharp frac- 
tional and full revivals. 
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FIG. 1: Schematic illustration of bistable wave packet mo- 
tion: the wave packet starts out in one diabatic state, splits 
at the crossing point and returns to the crossing where the 
parts interfere to reproduce the initial wave packet after one 
classical period of an individual potential well. The upper 
plot (a) displays the example of two bound states, typical 
of our trapped ion system, while the lower one (b) shows a 
bound state coupled to a repulsive one, which is a system of 
relevance to molecular dynamics (like (a)), where x represents 
the internuclear distance. 



Here we extend our former work and explore the dy- 
namics of a system describing a trapped ion driven by 
a standing wave. We point out that also in such a sys- 
tem bistable motion exists and investigate the short- and 
long-term evolution of the system under this interference 
condition. 

The outline of the article is as follows. We begin Sec. 
Ull by presenting the trapped ion model Hamiltonian, as 
it is generally stated in quantum optics, and discuss pa- 
rameter regimes and the corresponding approximations 
- rotating wave, Lamb-Dicke, weak and strong excita- 
tion and adiabatic. It is shown how the Hamiltonian re- 
laxes to an effective one mimicking the acclaimed Jaynes- 
Cummings (JC) model of quantum optics [37]. The re- 
lation between frequently used bases are considered in 
Sec. Ill Bl for which we define the corresponding poten- 
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tial curves. This enables us to obtain a deeper intuition 
about the dynamics and physical properties of the cou- 
pled system. The analogous molecular electronic states 
with their potential energy curves and couplings are also 
discussed, briefly. The following Sec. IIIII is dedicated 
to our numerical results. First we define the concept of 
bistable motion and then show how it enters into the 
trapped ion model. For short timescales, we show in 
Sec. IIII Al that this type of dynamics oscillates between 
an entangled (can be made maximally entangled) and a 
disentangled state. We analyze the nature of the wave 
packet evolution around the potential curve crossing in 
terms of the ionic inversion, the autocorrelation and the 
von Neumann entropy. For very long times (Sec. IIII Bp 
we observe distinct wave packet collapse, fractional and 
full revival structures. Finally, Sec. IIVI gives a short 
summary. 



II. THE MODEL SYSTEM 

This section introduces the model trapped-ion Hamil- 
tonian. The problem is first formulated in terms of 
operator algebra, commonly used in the quantum op- 
tics community, and then recast in the conjugate vari- 
able representation frequently applied in, for instance, 
molecular physics. Within the latter representation, 
the obtained Hamiltonian resembles that of an artifi- 
cial diatomic molecule with two coupled bound electronic 
states. There are several parameters determining the 
particular form of the potential curves and coupling, 
namely, effective ion-laser coupling A (Rabi frequency), 
trap frequency laser frequency lul, ion transition fre- 
quency ^, LD parameter (laser wavelength) 77, and laser 
phase shift <ft. For a large LD parameter, the curves 
exhibit several crossings, whereas for small r] one may 
encounter only a few, one or even no crossings in the 
parameter range of interest. We will here be interested 
exclusively in the situation of a single level crossing of 
the two potential curves. 



In order to gain a deeper understanding of the evo- 
lution we discuss the approximations regularly applied 
in trapped ion models. These often result in analytically 
solvable effective Hamiltonians. In the LD regime, rj <C 1, 
we may expand the cosine function to obtain 



H 



-cr z + A 



cos(</>) — 77 sin(</>) (a/ + a 



+o( v 2 )- 

In this regime, the effect of the phase shift <j) is clear. For 
(j) — kir, where k is an integer, the ionic vibrational states 
are not coupled by the external laser (carrier pumping), 
whereas if k = 1/2,3/2, ... the vibrational states do 
couple (sideband pumping) [38[. In the latter case, the 
Hamiltonian ([2]) coincides with the Rabi model [28] . Ne- 
glecting the virtual processes (i.e., aa~ and a^a+ corre- 
sponding to simultaneous deexcitation/excitation of the 
field and the ion) one derives the regular JC Hamiltonian 
in the RWA 



Hjc = ooa) a 



(3) 



where in this case g = A77 and a + = |2)(1| and cr~ = 
|1)(2|. The JC model is analytically solvable and the 
eigenvalues of its corresponding Hamiltonian are known 
as (up to an overall constant) 



n = 0,1,2, ... . 



Ef c (n)=um± ^ A 4 ^ )2 +<7 2 n, 

(4) 

As pointed out, the effective RWA Hamiltonian, in 
principle, is only justified for cos((/>) = 0, whereas in 
general situations more sophisticated RWA approaches 
must be considered. Turning to an interaction picture 
with respect to the free field and ion, Hamiltonian ([2|) 
becomes 

H T = UtHUjT 1 = 



\cos(<fi) (< 



a+e- iAt + a-e iAt ) 



A. Trapped ion in a standing wave 

In the formalism of phonon creation a) and annihila- 
tion a operators the Hamiltonian reads (we use atomic 
units so that h = 1) 



H 



\cos[q(a) + a) + (j)]cr x , (1) 



with the ion-laser detuning A = Q — ujl and the pr- 
oper ators act on the two ionic levels |1) and |2), a z = 
|2)(2| - |1)(1| and a x = |1)(2| + |2)(1|. The LD parame- 
ter is defined in terms of the ionic mass m and laser wave 
number k as 77 = k^/muu/2, and the relative position of 
the trap in the standing wave laser field is determined by 



-7/Asin(0) (tfa-e-^-^ + 



(j+ae^ 



-A)t\ 



-7/Asin(0) (a+6--e- i ^ +A ) t + . 

(5) 

where Ui = exp (— iuoa)at) exp (— i^a z t). This intro- 
duces three different time scales of the interaction, sug- 
gesting an application of the RWA. It is known, how- 
ever, that the validity of the RWA does not solely de- 
pend on the involved time scales, but also on the rel- 
ative amplitudes of the coupling parameters [HI, [29j]. 
In particular, only for ion- field couplings much smaller 
than the ion level separation and vibrational spacing - 
\r]sm((f)) <C ^,cj and Acos((/>) <C - is the RWA 

justified. We then note: Firstly, the Lamb-Dicke regime 
seems to favor application of the RWA, which, in fact, 
we have verified by numerical calculation. Secondly, the 
angle <fi is of importance for the validity of the RWA and 
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together with the three characteristic time scales it be- 
comes clear that a straight-forward application of the 
RWA in the basis presented above is non-trivial. The 
first observation will be of importance further on when 
we introduce the adiabatic approximation. The second 
is indeed an interesting issue that was resolved by Wu 
and Yang [17]. Their idea is that also the third term of 
the expanded Hamiltonian ([2]) is considered as a "free" 
part, and these first three terms define the interaction 
picture and the proper time scales. Within this frame, 
the slowly oscillating terms describe a JC type of interac- 
tion between the field (vibrations) and "dressed" atomic 
states [39j differing from the |1) and |2) states. 

Yet another, but related, approach is to use the expan- 
sion m m 



cos[?7(a^ + a) + (j)} 



k=0 



k\ 



[o*/fc(n + *) + /*(" + 



where 



fk(n) = £ 



(h — m+ l) m 
(fc + l) m m! 



fc! 



Ljjri 2 ) 
(n + l) fc 



(6) 
(7) 



with (n — m + l) m = h(h — 1) . . . (n — m + 1), h = a) a, 
€k — 1/2 for k = and = 1 otherwise and L^(r] 2 ) 
are the generalized Laguerre polynomials. Imposing suit- 
able resonance conditions (|A — u\ <C |A — fco;|, for 
k = 0,2,3,4, . . .) between the constituting frequencies, 
we may for proper amplitudes of the couplings employ 
a RWA to regain a JC like interaction in terms of "de- 
formed" oscillator operators A = afi(n + 1). The A and 
define generalized boson operators obeying modified 
boson algebras and they specify, for instance, non-linear 
Fock and coherent states being eigenstates of A^A and 
A, respectively. Note that, in some sense, the new ladder 
operators can be viewed as describing a non-linear oscil- 
lator with amplitude-dependent frequency. Additionally, 
it is known that the anharmonicity of an oscillator greatly 
affects physical phenomena like collapse and revivals [42[ . 
Related algebraic techniques, in particular shape invar i- 
ance and SUSY, can therefore be a useful tool for describ- 
ing molecular dynamics where oscillators are in general 
anharmonic [43[. We remark, however, that most of these 
models J43j employ the standard RWA as it is used in the 
JC model, while it is clear that this is in general less jus- 
tified than for the JC model, due to varying distance 
between energy eigenvalues, for instance. 

We pointed out above that the relative size between the 
effective laser-ion coupling A and the two-level spacing ft 
is crucial for the validity of the RWA. Also the relation 
between A and uo is important. In general, for moderate 
phonon numbers (n < 10), the strong-excitation regime 
[24j , A ^> lj, invalidates the RWA and the LD approxima- 
tions, while the opposite weak-excitation regime [22[ fa- 
vors them. It should be emphasized, though, that, even- 
tually, the applicability of the approximations depends 



on all parameters, including ft and 77, and their mutual 
relations. 

Now, rather than by the boson creation and annihila- 
tion operators as above, we represent the Hamiltonian in 
terms of conjugate variables 



V 2muo V ; ' 



(8) 



V 



obeying the regular canonical commutator relation. In 
this nomenclature we find 



H 



rriu 



2m 



A 
~2 

Acos(/cx + <p) 



\cos(kx + <p) 
A 



(9) 

where k is the external laser wave number and m the 
mass of the ion. This representation of the Hamiltonian 
serves as the starting point of our analysis; for a given 
initial state, its wave packet will evolve on two coupled 
harmonic oscillators according to eq. ([9]). We note that 
in ^-representation a Fock state of the vibrations, h\n) = 
n\n), or a coherent state, a\a) = a\a), read 



<p n (x) = (x\n) 
cp a (x) = (x\a) 



2 n v / / 7m! 
/mco\ 1 / 4 



1/2 



H n (^/muux)e 2 , 



V 7T 



e -iMa)] 2 e - 



Here, H n (x) is the n:th Hermite polynomial. 



(10) 



B. Base representations 

To gain a deeper intuition of the problem it is in 
order to discuss various bases that will be used in 
the following. One natural representation is that of 
bare basis states, which are the eigenstates of the free 
Hamiltonian part {|n, 1), |n, 2)}, where n represents the 
n:th Fock eigenstate of the harmonic oscillator. In 
^-representation the bare basis states are written as 



(x\n, 1) = <p n (x) 



, (x\n,2) = <p n (x) 



A general bare state is given by i/ji(x)\1) or ip2 ( x )\2) for 
some normalized wave function ipi (x) . From the form on 
which the Hamiltonian ([9]) is presented, it is clear that 
the external laser field enters in the off-diagonal terms 
and hence couples these states. The bare states can be 
used to define bare potential curves as the diagonal po- 
tential elements of the Hamiltonian in the ionic basis |1) 
and 1 2). Thus, the bare potentials are two centered os- 
cillators shifted in energy with either A/2 or —A/2. 

The eigenstates of the full Hamiltonian are called 
dressed basis states and denoted by |x+) and |x~). We 
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have introduced the superscript ± since, for a wide range 
of parameters, the eigenstates come in pairs, which is true 
in particular when the RWA or the LD approximation or 
both have been imposed, as seen in ([4]) for the JC model. 
In these regimes, not necessarily only for the JC model, 
"simple" analytic solutions to the problem are available, 
whereas in the general case this is not true. Note that 
there is no obvious corresponding set of dressed potential 
curves. 

The above two bases are part of the the conventional 
terminology of quantum optics. We now turn to a more 
convenient representation of coupled molecular electronic 
states. For this purpose, we employ the ^-representation 
and rotate the Hamiltonian ([9]) around the (a x -f <7 z )-axis 
by the unitary operator U\ = (a x + a z ) giving 



2m 



\cos(kx 
A 
2~ 



A 

-Xcos(kx ■ 



(ii) 

This transformation swaps the off-diagonal terms with 
the diagonal ones of the last term in the Hamiltonian. 
The form of Hi has many similarities with the Hamilto- 
nian describing a diatomic molecule in a diabatic repre- 
sentation. The off-diagonal elements in ([TT]) are in poten- 
tial form, i.e., they do not contain differential operators. 
Here, the couplings are x-independent, which however 
is also frequently assumed in models used in molecular 
physics, e.g., the Landau-Zener model Note that 

in our ion-trap system x represents the spatial center- 
of-mass position of the ion in the trap, whereas in the 
molecular counterpart x is the inter nuclear distance. The 
operator U\ transforms the bare internal ionic states as 



-) = U 1 \l) = 



V2 



(12) 



-1 



which defines the diabatic states as (x\ip)\±) = ip±(x)\±), 
and the corresponding diabatic potential curves as the di- 
agonal potential matrix terms of ([TT]) . Thus, contrary to 
the bare curves, these are centered oscillators modified 
by ±A cos(kx + <p) and depending on the parameter am- 
plitudes very different types of diabatic potentials may 
be obtained. For example, if the oscillating cosine func- 
tion dominates over the harmonic oscillator (in the spa- 
tial range of interest), the system is semi-periodic and is 
best understood from Bloch or Floquet theory (45[. In 
the opposite limit of a small A and k we may expand 
the cosine as in (j2j) to find the two centered oscillators 
both shifted and displaced according to the added terms 
±A [cos(0) — 77 sin ((/>)£]. In this paper, we are interested 
in a regime intermediate between the two. 



Given the Hamiltonian in the form of ([9]), the trans- 
formation 



U 2 = 



cos(#) — sin(0) 
sin(0) cos(#) 



with 



tan(20) 



2Xcos(kx + < 



(13) 



(14) 



diagonalizes the last term of the Hamiltonian. However, 
since U2 is x-dependent it does not commute with p and 
explicitly we get 



#2 = U2HU2 1 



e(x) 



rmo 



2m 2 

d 2 - 2i(d0)p 



{dey 



-d z 6 + 2i(d6)p 



where df = df/dx and 



-e(x) 



(15) 



e{x) = a/A 2 /4 + A 2 cos 2 (kx + </>). 



(16) 



These diagonal terms define the adiabatic potential 
curves, while U2 acting on the bare states gives the adi- 
abatic states; (x\i/j)\i) a d = i/j(x)\i) a d = VK^)^ K), for 
i = 1,2. Note that \i) a d is m general x-dependent. The 
diagonal and non-diagonal terms containing dO and d 2 
are the non- adiabatic corrections, and if these play a mi- 
nor role in the wave packet evolution, the dynamics is 
said to be adiabatic. For the present system, we have 



dO 



d 2 



AkXs 
"A 2 + A2c 2 ' 

A[Afe 2 c(A 2 +4A 2 c 2 )-8A 3 feW] 
(A 2 +4A 2 c 2 ) 2 



(17) 



where s = sin(kx + <p) and c = cos(kx + 0), and typically 
the second, higher derivative term, is smaller. For di- 
atomic molecules an equivalent transformation from the 
diabatic to an adiabatic representation can be applied. 
Consequently, also then non-adiabatic couplings are both 
x- and p-dependent. We remind, as discussed in the pre- 
vious section, that the RWA and LD approximations are 
highly related in terms of validity ranges and especially 
r\ <C 1 (or likewise for small k) is a constraint on their 
implementation. Here, as well, we note that for small 
k the non-adiabatic coupling terms (jT7j) become small. 
On top, we claimed that A < A is also a condition for 
the applicability of the RWA and the LD approximation, 
which from ([TT]) again is seen to favor the adiabaticity. 
We have verified by numerical calculation that all three 
approximations more or less share the same parameter 
validity range. From (fT7|) one notes that also the case 
when A = gives vanishing non-adiabatic corrections. 
This limit, however, does not represent a proper adia- 
batic one, but rather the diabatic limit, and is hence 
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State 


Definition 


Potentials 


Bare 


ipi(x)\i) 




Dressed 


Eigenstates of full Hamiltonian (|9]), \xt) 




Diabatic 


i/>-(x)\-) = tl>-(x)Ui\l), i/>+(*)\+) = i/>+(x)Ui\2) 


V£(x) = ^-x 1 ±Xcos(kx + 0) 


Adiabatic 


ipi(x)\i)ad = ipi(x)U 2 \i) 


V£{x) = ^x 2 ±s(x) 



TABLE I: States and potentials used in this article. Here, i = 1, 2 and ipi(x) is a general normalized wave function. 



sometimes referred to as the anti- adiabatic limit. Since 
we consider parameters where the adiabatic approxima- 
tion breaks down, it follows that we can apply neither 
the RWA nor the LD approximation. 

In Table I, we summarize the various bases, while in 
Fig. [2] we show some typical examples of the correspond- 
ing potentials. Note that at the crossing point, when the 
cosine function is zero, the bare and adiabatic potentials 
coincide, while away from the crossing the diabatic and 
adiabatic curves approach each other for the current set 
of parameters. We may point out that different defini- 



tions of the above bases exist, and in particular in quan- 
tum optics is the bare-diabatic and the dressed-adiabatic 
representations often identical. Identifying bare states 
as diabatic ones and dressed states with adiabatic ones 
does not, however, properly coinside with the common 
definitions in molecular and chemical physics. Adiabatic 
states are in general not eigenstates of the full Hamilto- 
nian, as dressed states are, and diabatic states usually 
approximate the adiabatic ones far from curve crossing 
regions, which is not the case of bare states. 




-6 -4 -2 2 4 6 
X (a.u.) 



FIG. 2: Example potentials: bare (crosses), diabatic (solid 
line) and adiabatic (dots). In (a), the potentials only possess 
one crossing in the shown interval, while (b) displays an ex- 
ample of several crossings. The parameters (in atomic units) 
are: m = 80000, uj = 0.0005, A = 0.05, A = 0.02514 in both 
plots, while k = 0.2 and = 1.07249074 rad in (a) k = 1 and 
= rad in (b). 



III. DYNAMICS - BISTABLE MOTION 

The existence of bistable evolution is a direct con- 
sequence of interference between the constitutent wave 
packet parts, as shown in Fig. [TJ The bistable motion 
turns out to be stable over very long time periods, show- 
ing clear oscillations, collapses and revivals [36[. In or- 
der for the wave packet fractions to overlap and inter- 
fere correctly, a necessary, but not sufficient, condition 
for bistable motion is that the two parts of the splited 
wave packet should return simultaneously to the starting 
point. We denote this period of time by T c ^, where the 
index cl stands for classical, and it is determined from the 
first maximum of the autocorrelation function, as will be 
described below. Further, let P%(t) be the population in 
diabatic state \i) (we hence trace out the field degrees of 
freedom). If, initially, the wave packet is in the diabatic 
state i we define P sp = 1 — Pi{t = T c //2), which gives 
the population in the opposite diabatic state at half a 
classical period. In other words, P sp in a sense measures 
the amount of splitting of the wave packet as it traverses 
the crossing and, hence, for times shorter than the col- 
lapse time Pi(jT c i) « 1 (j = 0,1,2, . . .) for a bistable 
trajectory. Naturally, < P sp < 1, where P sp = cor- 
responds to diabatic and P sp — 1 to adiabatic evolution, 
and in particular, bistable trajectories are not restricted 
to a certain P sp . It is clear, however, that for the non- 
trivial intermediate values of P spi these types of inter- 
ferences are sensitive to parameter values. In addition 
to bistable trajectories, the parameters may be tuned so 
that P»((2j + l)T d ) w and Pi(2jT d ) w 1, where again 
j = 0, 1, 2, . . ., which defines an astable trajectory. Here, 
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however, we focus on the more stable bistable case. 

For a general set of bound coupled (or uncoupled) po- 
tential curves we may expand the eigenenergies around 
the mean set of quantum numbers no forming the wave 
packet. For a single discrete number n this reads 



E(n) « £(no)+£'(no)(n-n ) + 



E"(no) 



(n - n ) 2 



E'»(no) 



(n - n ) 3 + 



(18) 

where E'(no) = (dE(n) /dn) n=no and so forth, and no ^> 
1 is the "average" quantum number of the wave packet. 
The wave packet is assumed to be composed of eigen- 
states with quantum numbers n fairly localized around 
no- We define the different time scales 



2tt 



2tt 



2tt 



|£'(n )|' 



• J ~ |£"(n )|/2' ■'-sup - |£y"( no )|/6' 

(19) 

characterizing the classical vibration, revival and superre- 
vival time, respectively. Note that for a single harmonic 
oscillator T c \ = 2tt/u; and the higher order terms van- 
ish, which is equivalent to letting the latter characteris- 
tic time scales go to infinity. For two uncoupled different 
harmonic oscillators 1 and 2 one has the respective clas- 



sical periods 



(2) 

2tt/uji and = 2tt/uj2, where uji 



is the frequency of the oscillator i. The combined clas- 



sical period becomes T c \ — T^k = T^U, where k and I 
are the smallest possible integers obeying the condition 
kuo2 = Thus, it is clear that for multi- level systems 
the typical time scales can be very long. As an example, 
returning to the JC model (|4j), the revival time is usually 
defined in the quantum optics community as the time 
it takes for the constituent neighboring Rabi frequencies 
Q n = v/(A - c^) 2 /4 + # 2 n to differ by tt, 



-.(2) 



(a 



no + l 



^no) T rev — TT. 



(20) 



Consequently, the JC revival time as it is given in ([20]) 
is more reminiscent of the classical period rather than 
the revival one according to the definition ([T9|) . Also the 
fractional and superrevivals as discussed in the literature 
of the JC model differ from the ones of eq. ([T9]h see [46| . 

The wave packet propagation is carried out using the 
Chebychev polynomial method [47[. We assume the ini- 
tial state to be of the form ^(x,0)|+), where ip(x,0) is 
taken to be a minimum uncertainty Gaussian function 



*M)=#z,o)|+) 



27RT 2 



1/4 



(x-x ) 



(21) 



T rev . The quantities of interest for us are 
W{t) = ((T z ), Inversion 



A(t) 



V*(x, t)#(x, 0)dx, Autocorrelation (22) 



S T (t) = -Tt [ Pi log(p/)] , Entropy 

where ^(x,t) is the full two-level wave packet at time t, 
and pi = Tip [p] is the reduced density operator for the 
ion, obtained once the field degrees of freedom of the full 
density operator p = \^f(x,i))(^f(x,i)\ have been traced 
out. Note that W(t) is the inversion between the bare 
states 1 2) and |1). The autocorrelation function deter- 
mines the overlap between the initial state with the one 
at time t, while the von Neumann entropy is a measure of 
the degree of entanglement shared between the field and 
the ion. For initial pure states, as in our case, the entropy 
is the same for the field and the ion; Si(t) = Sir(£) |48j |. 
An advantage of our numerical wave packet approach is 
that, once the full wave packet ty(x,t) is given, all the 
above quantities are easily calculated, as well as other 
relevant properties. 



A. Short time evolution 

By short times we refer to time scales of order T c \ <C 
T rev corresponding to a single or few oscillations of the 
total wave packet. We tune the parameters to have 50- 
50 splitting of the wave packet at the curve crossing, i.e. 
P sp = 1/2, corresponding to maximal entanglement. The 
initial wave packet is localized around a = xq yj moo / \/2 = 
Qyjmuj j \f2 ~ 26.8 for the considered mass, and thus for 
the chosen parameters the motion is highly excited with 
a mean phonon vibrational number no = 720. Its width 
is either a = 0.047 (highly squeezed |50[), or we pick a as 
for a coherent state with given mass m. In atomic units 
m = 80000, which corresponds to 42 amu. The results 
are applicable to any mass, however, provided the time 
is scaled accordingly. The rest of the parameters are as 
in fig. [2] (a), except A that should be divided by 5. 

In fig. [3] we display the numerical results of the quan- 
tities ([22]) , (a)-(c) refer to an initial coherent state and 
(d)-(f) an initial Gaussian distribution with a = 0.047. 
We note that in both cases almost perfect 50-50 splitting 
is obtained and thus the system is maximally entangled 
on one side of the level crossing, while it is disentangled 
on the opposite side. The state can be written as 



= tl>+(x,t)\+)+il>-(x,t)\-) 



V2 



OMM+^+0M))|i> 



(23) 



including also the case of coherent states. This wave 
packet is let to evolve over long time periods, typically 



+ _(^_(z,t)-<MM)|2). 
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For 50-50 splitting, and well separated splitted wave 
packets, we have that the vibrational wave functions 
of the bare states |1) and |2) are orthogonal, in agree- 
ment with the maximal entanglement. At this instant, 
the combined bi-partit ion-phonon system is in an EPR 
Schrodinger cat state. To a first approximation, we may 
assume the wave packets to keep their shape over a time 
period of one classical oscillation. That is, for initial co- 
herent states, the splitted states are still coherent with 
modulated phase and amplitude. This is indeed verified 
below, where we considers the dynamics in phase space. 
Thus, the phonon distrubutions in this case are Poisoni- 
ans for the diabatic states, while for the bare states they 
would be made out of the sum of two coherent states with 
equal phase but different amplitude. 

Calculating the variance of the wave packets, one notes 
that the squeezed state [50|, a = 0.047, shows typical 
breathing motion while traversing the potential back and 
forth. From the decrease in amplitudes, the quantities 
corresponding to the coherent state is seen to be more 
sensitive than for the squeezed state, or in other words, 
the collapse time is faster. This fast dephasing comes 
about due to the larger spread in phonon numbers n for 
the coherent initial state. It also shows up as rapid os- 
cillations in the inversion around the curve crossings. 

The plots also give a measure of the characteristic tran- 
sition time at the level crossing [5lj, and one notes that 
this time is rather short compared to the full time X^/, 
which is indeed known from other similar models [521 ]. 
Interestingly, we have found that short crossing time as 
well as the sharp peaks with large amplitudes in the 
inversions, (a) and (d), around the crossings are typi- 
cal for the bistable cases. The bare state inversion can 
be expressed in terms of the diabatic wave packets as 
(a z ) — 23? [J ip^_(x 1 t)ip-(x 1 t)dx] 1 and is hence related to 
the coherence of diabatic wave packets. We note that the 
total energy can be written as 

E to t{t)= [^Ux,t)Hu^x,t)dx + jW(t), (24) 

*=+,- 

where and H are the two diagonal matrix ele- 

ments of the Hamiltonian ([TT]) . From this it is seen that 
the bare state inversion is a measure of the diabatic en- 
ergy transfer. 

At first, it may seem surprising that the inversion W(t) 
can be non-zero at the right side of the crossing, even if 
the autocorrelation A(t) is close to 1 and the entropy 
Si(t) is close to zero. The zero entropy, however, is not a 
contradiction as the field and ion states can be separated, 
provided the constituent wave packets fractions ijj+(x : i) 
and ip-(x,t) differ merely by a complex constant. For 
the inversion we have 



w(t) = (Mt)\Mt)) - (Mt)\Mt)) 

= 2R[(V>+(f)|V-(f)>] 



(25) 



while for the autocorrelation 

\A(t)\* = |(*W|*(0))| 2 = |^ + W|^ + (0))| 2 

< l-<V»_(t)^_(t)>. 



(26) 



For sufficiently small (il;-(t)\ip-(t)) it follows that A(t) 
can be close to one even if W(t) is not exactly zero. 
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FIG. 3: Inversion W(t), autocorrelation A(t) and entropy 
Si(t) for a coherent (a)-(c) and a squeezed (d)-(f) initial state. 
See text for the parameters. 



The wave packet splitting manifests itself also in the 
phase space distributions [281 ]. In Fig. [4] we display snap- 
shots of the Wigner distribution [53[ 



(2tt)2 



df / dr e^+^TYj 



e -i(rp+£x) 



- dye 

7T ' 



-i2py 



-H/>* (x-y)i/;-(x + y) 

(27) 

over one period of oscillation for a coherent state (a) and 
a zoom in at the Wigner distribution close to half an 
oscillation for a coherent (b) and squeezed state (c). In 
(d) we show the Wigner distribution after the collapse 
at approximately t « 2000T C ^ for a squeezed state. The 
spreading of the distribution over the whole accessible 
phase-space is a character of the collapse. At times of 
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(a) 





fractional revivals, the wave packet forms multiple local- 
ized bumps [5J]. The two constitute parts of x) 
approximately follow the "classical" phase space trajec- 
tories 



dxc 
dt 



dp c 



dp c 
dt 



8Xr 



(28) 



with the adiabatic and diabatic Hamiltonians = 
£s_ _|_ yiv ( x ) and i = A, D. The small anharmonicity 
of the potentials is reflected in the pattern of the Wigner 
distribution. The squeezed distributions rotate along the 
trajectories, which gives rise to the characteristic breath- 
ing motion of its width. The phase space plots suggest 
that the two involved wave packets roughly evolve on 
either of the potential curves VP{x) or V+(x), and at 
the right side away from the crossing point these two ap- 
proximately coincide. We argued that the corresponding 



classical times must be equal, T cl + 



in order for 



the wave packets to overlap at the curve crossing and in- 
terfere maximally. This, however, does not imply that 
the longer time scales, T^ ev , T® ev1 T^ p , . . ., of the poten- 
tials are equal, as will be discussed next. 



-200 



,03, 
CL 



200 




B. Long time evolution 



-200 



We now turn to the dynamics over time scales of the 
order of T revi which for our parameters correspond to 
several thousands of classical oscillations. Since we are in 
the strong-excitation regime, this does not, however, nec- 
essarily imply exceedingly long total operational times. 

The expansion of the eigenvalues ([T8]h in general, con- 
tains an infinite number of terms, but for relatively 
smooth potentials one expects fairly fast convergence. 
Higher order terms typically cause imperfect full time 
and fractional revivals, by which we mean that the am- 
plitude of, say, the autocorrelation function squared does 
not reach the values where j = 1,2,3, . . ., as 

expected [Hj for rational fractions of the revival time 
kT rev /j, k and j mutually prime. This effect, arising 
from higher order anharmonicities, is naturally more pro- 
nounced for longer times, and hence the revival ampli- 
tudes drop for each consecutive revival period T rev . The 
widths of the revival peaks are getting broader, as well. 
As seen in Fig. [3j the uncertainty in phonon numbers n 
affects the collapse-revival pattern. In particular, the en- 
velope functions of the revival peaks are highly sensitive 
to the uncertainty An, while the location of the peaks are 
determined by the average no- In the "ideal" case, full 
amplitude revivals occur at t = jT rev /2 for integer j, half 
amplitude revivals at t = (2j + l)T rev /4, and so forth. 
In Fig. [5]we present the autocorrelation function (a) and 
von Neumann entropy (b) for a time span slightly longer 
than T rev . The parameters and initial conditions are sim- 
ilar to the ones of Figs. [3(d)- (f); = 1.07244080531656 
rad, A = 0.064727653164347, a = 0.0340999659 and 
A = 0.005025343787836. The slight change in param- 



2 
Oh 



200 



-6 



-2 2 
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FIG. 4: (Colour online) The Wigner distribution for an initial 
coherent state (a) and (b) and for a squeezed state (c) and 
(d). The first plot shows snap-shots throughout one period 
of oscillation, while (b) and (c) zooms in at the distributions 
around half oscillation time. In (d) the distribution at the 
time of collapse is displayed. 



eters gives 40-60 splitting (P sp = 0.4) rather than 50- 
50. Some of the numerous fractional revivals are labeled 
conventionally [42[. From Fig. [5j(b), we note that in the 
collapse region the field and ion are highly entangled, but 
the fact that P sp =0.4 causes a decrease of entanglement 
compared to the P sp = 1/2 case. 

We have previously [36[ studied the bistable motion 
and long-term revivals in system Hamiltonians with more 
resemblance to molecular models. In terms of the clas- 
sical period, the revival time found in the present study 
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is much longer than those we calculated for the molecu- 
lar systems. This originates in a smaller anharmonicity 
of the ion-trap system. In [36[, we found that in accor- 
dance with expectations from a semiclassical derivation 
the revival time of the coupled system T rev can be ob- 
tained from the revival times of the adiabatic T^i and 
diabatic T^v pathways and the wave packet splitting as 



1 



1 



-I- rev l rev 



rpD- 



(29) 



This relation was fulfilled for numerous systems, com- 
prising both coupled bound-bound and bound-repulsive 
states. Equation ([29]) accurately reproduces the limit- 
ing cases of adiabatic P sp = 1 and the diabatic P sp = 

evolution. If T^t ^ T^ev , then we presumably have 
that the size of T rev is between the values of T re v and 
T^ e v in the intermediate range. However, in reality T^t , 
Tr° e y and P sp are not independent variables as they all 
depend on the system parameters in a complex manner. 
To conclude, the revival time for multi- level systems can 
be shorter than those of individual isolated systems. For 
the classical periods, on the other hand, we found that 
T c i = = for bistable trajectories. Thus, the 

dynamics cannot be viewed as two uncoupled wave pack- 
ets evolving on the potential curves V+{x) and VP(x). 
Only interference between these two wave packets can 
cause a common revival time T rev shorter than one of 

Trev Or Trev • 
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FIG. 5: Autocorrelation function (a) and von Neumann en- 
tropy (b) for an initial squeezed state in the case of 40-60 
splitting. 



IV. CONCLUSION 



In this article we have studied in a non-standard way 
the dynamics of a harmonically trapped ion pumped by 
a standing wave beyond the RWA and the LD approxi- 
mation. Using a wave packet technique, we predict the 
existence of bistable states of motion similar to those ob- 
served for molecular systems [33|, [34], HH, [36| . These arise 
from interference of wave packets and hence it is a pure 
quantum effect. Over longer time periods these states 
possess a well resolved collapse-revival pattern including 
full and fractional revivals. 

Beyond the conjecture of the new class of bistable ion- 
trap states, the article establishes a link between trapped 
ion systems and molecular physics models. As such, we 
have put extra stress on how the two areas relate to 
one another; expressing the model Hamiltonian in x- 
representation and introducing diabatic and adiabatic 
states and potential curves. We intend to further in- 
vestigate these directions, both by considering algebraic 
methods in simple molecular models and by extending 
the current ion-trap system to include additional ionic 
internal levels and vibrational degrees of freedom. We 
will also allow for different trapping potentials in order 
to more closely mimic molecular-like situations. Anhar- 
monic traps have been discussed in the literature, and 
then mainly quadrupole trapping [55], but in principal 
more exotic shapes may be gained by reforming and ad- 
justing the geometry of the trap electrodes. On top, in- 
dividual trapping potentials for the internal ionic levels 
have been discussed [7[. 
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